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Abstract — The underwater acoustic channel is characterized 
by a path loss that depends not only on the transmission 
distance, but also on the signal frequency. Signals transmitted 
from one user to another over a distance I are subject to a power 
loss of l~ a a(f)~ l . Although a terrestrial radio channel can be 
modeled similarly, the underwater acoustic channel has different 
characteristics. The spreading factor a, related to the geometry of 
propagation, has values in the range 1 < a < 2. The absorption 
coefficient a(f) is a rapidly increasing function of frequency: it 
is three orders of magnitude greater at 100 kHz than at a few 
Hz. Existing results for capacity of wireless networks correspond 
to scenarios for which a(f) = 1, or a constant greater than 
one, and a > 2. These results cannot be applied to underwater 
acoustic networks in which the attenuation varies over the 
system bandwidth. We use a water-filling argument to assess the 
minimum transmission power and optimum transmission band 
as functions of the link distance and desired data rate, and study 
the capacity scaling laws under this model. 

I. Introduction 

The seminal work by [1] studied wireless networks, mod- 
eled as a set of n nodes that exchange information, with the 
aim of determining what amount of information the source 
nodes can send to the destination as the number n grows. 
The original results obtained for nodes deployed in a disk 
of unit area motivated the study of capacity scaling laws in 
different scenarios, ranging from achievability results in ran- 
dom deployments using percolation theory [4] or cooperation 
between nodes [3], to the impact of node mobility over the 
capacity of the network, e.g. [2]. Reference [5] provides a 
good overview of the different assumptions and scaling laws 
for radio wireless networks. 

The underwater acoustic channel is characterized by a path 
loss that depends not only on the transmission distance, but 
also on the signal frequency [7]. Signals transmitted over a 
distance / are subject to a power loss of l~ a a(f)~ l . Although 
a terrestrial radio channel can be modeled similarly, the 
underwater acoustic channel has different characteristics. The 
spreading factor a, related to the geometry of propagation, 
has values in the range 1 < a < 2, where a = 1 corresponds 
to cylindrical spreading. Also, the absorption coefficient o(/) 
is a rapidly increasing function of frequency, e.g. it is three 
orders of magnitude greater at 100 kHz than at a few Hz [7]. 
Finally, the power spectral density of the noise underwater is 
highly dependent on frequency. 

Existing capacity scaling laws for wireless radio networks 
correspond to scenarios for which a(f) = 1, or a constant 



greater than one, and a > 2, e.g. [1], [4]. These results cannot 
be directly applied to underwater acoustic networks in which 
the attenuation varies over the system bandwidth and a < 
2. We study the scaling laws under a model that considers 
a water-filling argument to assess the minimum transmission 
power and optimum transmission band as functions of the link 
distance and desired data rate [8]. In particular, we study the 
case of arbitrarily deployed networks in a disk of unit area, 
and follow a similar procedure as in [1] to derive an upper 
bound on capacity. In this sense, we provide an extension of 
the work in [1] under a more complicated power loss model. 

We show that the amount of information that can be ex- 
changed by each source-destination pair in underwater acous- 
tic networks goes to zero as the number of nodes n goes to 
infinity. This occurs at least at a rate n^ l / a e^ w ^ 0< ^ n 
where Wq represents the branch zero of the Lambert func- 
tion. We illustrate that this throughput per source-destination 
pair has two different regions. For small n, the throughput 
decreases very slowly as n increases. For large n, it decreases 
almost as -nT l l a . Thus for large enough n, the throughput 
decreases more rapidly in underwater networks than in typical 
radio channels, because of the difference in the path loss 
exponent a. 

The paper is organized as follows. In Section II, we present 
the underwater channel model. In Section III, we analyze 
the scaling laws for the case of a network transmitting in 
an arbitrarily chosen narrow band. In Section IV, we study 
scaling laws for the low-power/narrow-band case, with optimal 
bandwidth allocation using a waterfilling argument. In section 
V, we consider the case in which the nodes can transmit at 
high power over a wide transmission band. Conclusions are 
summarized in Section VI. 

II. Underwater Channel Model 

An underwater acoustic channel is characterized by an 
attenuation that depends on the distance / and the signal 
frequency / as 

A(lJ)=(-^-Ya(f) 1 (1) 

V l ref J 

where l re j is a reference distance (typically 1 m). 

A common empirical model used for the absorption a(f) 
is Thorp's formula [7] which captures the dependence on the 
frequency. This absorption a(f) is an increasing function of 
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Fig. 1. Relationship between transmission distance and center frequency in 
a narrow band system. 



/. The spreading factor describes the geometry of propagation 
and is typically 1 < a < 2, e.g. a = 1 and a = 2 correspond 
to cylindrical and spherical spreading, respectively. The noise 
in an acoustic channel can be modeled through four basic 
sources: turbulence, thermal noise, shipping, and waves. It has 
a power spectral density (psd) which depends on the frequency, 
the shipping activity s, and the wind speed w in m/s [7]. 

The complete model for a colored Gaussian underwater link 
was presented in [8] where power was allocated through wa- 
terfilling. In the absence of multipath and channel fading, the 
relationship among capacity, transmission power, and optimal 
transmission band of a point-to-point link is given by [8] 



C= / log 2 

<B(l,C) 



K(l,C) 
A{lJ)N(f) 



df 



(2) 



where N(f) is the psd of the noise, B(l, C) is the optimum 
band of operation and K (I, C) is a constant. The transmission 
power associated with a particular choice of (I, C) is given by 



P(l,C) 



B(l,C) 



S(l,CJ)df 



(3) 



where the psd of the signal is S(l, C, f) — K(l, C) ~ 
A(lJ)N(f),feB(l,C). 

A distinguishing feature of the underwater acoustic channel 
is the dependence of the optimal transmission band on the link 
distance [8]. Fig. 1 illustrates the optimal center frequency 
f c (l) as a function of distance. The optimal center frequency 
is defined as the frequency at which A(l, f)N(f) is minimal. 
This implies that if the transmission power for a link is 
low, the transmission bandwidth will be low and around the 
optimal frequency. Thus, the optimal transmission band in the 
spectrum changes dramatically with the link distance. Fig. [T] 
also illustrates that a node transmitting over a short range 
will optimally be assigned a transmission band at high center 
frequency, as in case (a), while a node transmitting over a 
longer distance will be assigned a different transmission band 
at lower center frequency, as in case (b). 

For the case in which power available for transmission is 
low, the bandwidth of the transmission band will also be small. 



When the bandwidth is low enough \B(l, C) \ = Af, such that 
the product A(l, f)N(f) does not change much over that band, 
one can make a Taylor series approximation around the center 
frequency f c (l). This allows us to determine the power P for 
which the transmission band is narrow due to our waterfilling 
argument. The Taylor series approximation has the form 



A(l,f)N(f)^A(lJ c )N(f c ) + T 



(/ - fcY 



(4) 



V/ in the band, where T = fp (A(l, f)N(f)) | /=/c . Substi- 
tuting the expression Q into expression ([3]), and using the fact 
that K(l,C) = A(l,f mas .)N(f max ) = A(lJ min )N(f min ), 
where / max and / m j n are the maximum and minimum fre- 
quencies of the transmission band, we obtain 



P(Z,C)« A(Z,/ max )7V(/ max )A/ 



it:U(Uc)N(f c)+ 



•U-fcY 



df 



(5) 



where Af = / max - /min- Considering / max - f c « 
and f c — /min ~ -J- given our quadratic Taylor series 
approximation of A(l, f)N(f), the above expression reduces 
to 



T 

P = —A/ 3 
12 J 

III. Fixed Narrow-band model 



(6) 



Let us study the physical model [1] to obtain an upper bound 
to transport capacity for transmissions in a arbitrarily chosen 
narrow band in an underwater channel. The narrow band 
assumption allows us to consider the attenuation as a constant 
over that band. We assume that the nodes are arbitrarily 
deployed in a disk of unit area as in [1], that each node has 
an intended destination node, and that the requirement for 
successful reception at node j of a transmission from node 
i is 



A(\Xi-X m \,f) 



N(f) + EfceT,*#< A{\x^-x m \j) 



>/3 



(7) 



where X\\ is the position of node i, is the position of 

node j to which i is transmitting, and t is the set of all nodes 
transmitting simultaneously in the same transmission sub-band 
and time slot. We assume that all sub-bands are in the narrow 
band [1], so that the attenuation is only dependent on the 
central frequency of the narrow band. The parameter / is kept 
to keep in mind the frequency dependence, and to allow us to 
use these results in the following sections where we analyze 
more complex settings. We consider that A is the throughput 
of each node, the network transports XnT bits over T seconds, 
and that the average distance between source and destination 
of a bit is L. Finally, we define the transport capacity as XnL 
bits-meters per second [1]. 

Let us define W — A/log2(l + /3) to be the transmission 
rate, where A/ is the bandwidth of the narrow band chosen 



for transmission. Since \X). — Xj^\ < -4=, and a(/) > 1,V/, 
then the path loss is 



A(\X k -X m \,f) 



< 



ref 



a(f) 



2/0F = 



Using a similar procedure as in [1], we have that 

A(\X X I n< + 1 7A(/) 



(8) 



(9) 



Let us sum over all transmitters { £ r and use the definition 
of the path loss in expression ([TJ: 



^|^-^ (4) ra(/)^- X ^)l< 7( 



fl- 



ier 



fl 



where 7 a = '" e j7 = a (/) 2 ^ v/ ^- Summing over all 

sub-bands and time slots and dividing both sides by H, we 
obtain 



-^^r (h,b) < la — — 



H 



=1 h=l 



(11) 



where represents the h-th hop of a bit 6, and if is 

defined as the number of hops performed in T seconds, which 
can be bounded by H < 2^ [1]. Finally, r A (h,b) = 
l(h,b) a a(f) l ( h ' b \ where l(h,b) represents the distance be- 
tween receiver and transmitter for the h-th hop of bit b. Since 
the function r A (l) = l a a(f) 1 is increasing and convex for 
I > 0, a > 1 and o(/) > 1, then 

' ,.s\nTh(b) \ a I \nTh{b) 



6=1 h=l 



6=1 h=l 



<(lna(/))% 



3+1WT 
'~f3 iT 



(12) 



Let us define i/j = (In a(f)) a j a ^ ^f, and note that ^ > 0. 
The left hand side of the above inequality is a Lambert 
function of the form W a exp W, which is an increasing 
function when W > 0. Then, 



lna(/) 
H 



XnTh(b) I / ,l/ a \\ 

53^i(^6)<^/«exp^-Wb(^-JJ (13) 



where Wo(-) is the branch zero of the Lambert function, using 
the nomenclature in [6]. This fact implies that 



Substituting for ip in (14) , we obtain 



(14) 



q-l 

Since ff ~a is an increasing function for a > 1, and constant 

^ a— 1 

for a = 1, then < {^^j a . Considering that 

Wq(-) is an increasing function, we have that 



W 



\ja" 



> W 



2\na(f)a(f) a ^ ((3 + 1 \ 1/a 2 1 / a 



A/a 



Using these inequalities into expression ( |15) , we obtain the 
scaling law: 



XnL < QWn'^r exp ( -W ( $ 21n ^^ Q X 



(16) 



(10) where 



2 \ 1/a 

o(/)V5? 



Since the zero-branch of the Lambert function satis- 
fies Wo(x) > 0,Vx > 0, then the exponential term 
exp (— Wq {o ( n nas va l ues between and 1. Note 
that as n — > oo, the exponential term in the scaling law goes 
to 1. This implies that the exponential term is important to 
determine the scaling for n small, while for large enough n 

a — l 

the upper bound is 0{n a ). 

If we consider a(f) = 1, i.e. the same path loss model as 
in [1], and recall that Wq(0) = 0, then 



XnL < 



1 (2(3 + 2 



l/a 



Wn 



(17) 



- ( 8 + 1 \ 1/a ( (^ a W 

XnL < (^—W) ex P l-W M (15) 



which is the original result of [1]. We have thus proved that 
the result in [1] is valid for a > 1. 

Fig. [^illustrates the upper bound on XL for different values 
of a(f) ranging from 1 to 10,000, which are characteristic 
of an underwater environment at different frequencies. For 
example, a(f) = 1,000 corresponds to a frequency of around 
100 KHz. We have used a—l and the parameters specified 
in the figure. We also plot dashed lines proportional to n~ l l a . 
As expected, as n is large enough the exponential term of the 
upper bound becomes negligible making the bound scale as 
0(n -1 / ). However, for small values of n the bound starts at 
a common point for the different a(f) values, and decays very 
slowly. Fig. [5] also illustrates that the value of a(f) determines 
the transition between these two operating regions: the larger 
a(f), the greater n has to be before transitioning. Of course, if 
we use a transmission band with high a(f) each node will have 
to be able to transmit at higher power to reach its destination. 
In the underwater channel, it also means that a higher center 
frequency is required because a(f) is an increasing function. 

Finally, Fig. [2] shows that XL remains almost constant for 
n < 100 nodes, a(f) > 100 and a disk area of 1 km 2 , which 
corresponds to densities of up to 100 nodes per km 2 . Note 
that the expected density of nodes in an underwater network 
is usually much lower given the applications for which they 
are deployed, e.g. environmental measurements. Thus, XL is 
almost constant for practical purposes. 
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Fig. 2. Upper Bound on XL for an arbitrarily chosen narrow band and 
different values of a( f). Parameters used are W = 1 bps, a = 1, /3 = 2, 
area = 1 km 2 . 



For each of the different bands, the analysis is as before up 
to equation ( fT0| ). At this point, let us define ja(frn) as j a 
for band m. Let us sum over all sub-bands and time slots to 
obtain 



X j(i)\ a (fm) 



< 







-AWT J2 7« Urn) 



(18) 



mer 



where 5 is the set of time slots. We can use the fact that 
a(/m) > a min , where a min = min mgr a(f m ). In the uncer- 
water scenario, a m j n = a(/ m j n ) because a(/) is an increasing 
function. Defining r A (h,bJ min ) = l{h,b) a a{f min ) l ^ b \ 
and following similar steps that lead to Eq. (jTTJ we get 



\nTh{b) 



H 



}^}^r A {h,bJ min ) < 



(3 + 1 AWT 



6=1 h=l 







mer 



la(fm) 



Defining V = (lna(/ min )) a ^^ E me r7a(/m), we 
can use a similar procedure as in the previous section, to show 
that the scaling law now becomes 



IV. Low power - Narrow Band Case 

As mentioned in the introduction, one of the characteristics 
of the underwater acoustic channel is that the optimal transmis- 
sion band using the waterfilling principle depends strongly on 
the distance of a link. In particular, if the transmission power 
of a node is very low, then nodes will optimally transmit in dif- 
ferent bands corresponding to different transmission distances. 
Thus, interference will come only from nodes transmitting 
in the same band. We have derived an expression for the 
power under these assumptions in Section EL In order to assign 
disjoint transmission bands, we divide the total transmission 
band of the system into non-overlapping bands of bandwidth 
A/. We use f c (l) as the mapping between the transmission 
distance and the corresponding transmission band for a low- 
power/narrow-band scenario. Thus, if a node transmits to 
another node at a distance I, we assign the transmission band 
that contains the frequency f c (l) as in Fig. [TJ 

The capacity analysis is similar as before if a(f) is replaced 
by a(fm) f° r ea ch of the bands, where f m is the central 
frequency of transmission band m. Let us assume that each 
node is capable of transmitting at AW bps in each band, 
where AW = A/log 2 (l + 0), and Af is the bandwidth of 
each non-overlapping band. 

A. Multi-Node Hopping 

Note that the definition of H changes slightly when we 
allow multi-node hopping. In this case, H < T]F\AWn _ 
T ^ n , where V is the set of sub-bands used by the network, 
wdW=\T\AW. 



a-i ( ( 21na(f Trli J ( \^ 
XnL < <&Wn-^r exp -W $ ^ I * 



where 



mer / 



The scaling law is similar in structure to the one obtained in 
Section EH. However, the constant $ depends on the average 
of a function of the absorption coefficients at / m ,Vm € T 
instead of a particular value. Again, if a(f) = 1,V/ the result 
reduces to that of [1]. 

B. Direct Transmissions 

If we constrain our system to perform direct transmissions 
only, using the fact that there is an assignment of frequency 
bands in terms of the distance, we can consider that h(b) = 
1,V6, i.e. only one hop. Given the distance-band separation, 
the problem can be thought of as solving for several networks 
that lie on top of each other, in different layers with no cross- 
layer interference. Membership to the layers is based on the 
distance of the connection. In other words, each transmission 
band m will have n m transmitters, where n — X^mer n? " 
constitutes the total number of nodes in the network since 
each transmitter has only one intended destination. 

This causes a different capacity scaling for each of the 
transmission bands, i.e. the scaling for each transmission band 
will have the form of expression ( [To} with AW instead of W 
and n m instead of n to obtain the scaling for band m. 



V. High Power - Wide Band Case 

In this scenario, nodes have enough power to transmit in a 
wide transmission band B, which implies that the absorption 
cannot be considered to be a constant over the band. The band 
B is again chosen using a waterfilling argument. We consider 
that the SINR requirement can depend on the frequency. That 
is 



Pj(f) 
A(\Xi-X m \,f) 



P k U) 



>/?(/) 



(19) 



N (f) + Efc 6 r,fc/i A{\X k -X m \J) 

We define W as the transmission data rate over the entire 
band, computed as 



W = 



log 2 (l + /?(/))# 



(20) 



If we assign a transmission rate to every sub-band df of dW = 
log 2 (l + (3(f))df, the analysis for each frequency is similar 
as in Section IV by letting A/ — * 0, renaming A/ as df and 
replacing the sums by integrals. Then, we have that 



\nTh{b) 
6=1 h=l 



< 



T 
H 



(3(f) + 1 
(3(f) 



la{f)dW 



T 
H 



(3(f) + 1 



7a(/)log 2 (l + /3(/))d/ (21) 



IB P(f) 

where / m j n = argminj a(f) and H can be shown to have the 
bound H < T ^ n using the definition of W in expression |20|. 
Following the procedure of Section IV, we show that the 
scaling law for the high power - wide band case has the form 



XnL < QWn^r exp -W 6 



21no(/ n 



where 



e 



2 1/a ( 1 f ((3(f) + 1) a(f) h io g2 ( i + /?(/)) ^ 1 



B 



W) 



VI. Conclusion 



This work presents upper bounds on the transport capacity 
of underwater acoustic networks with nodes deployed arbi- 
trarily in a unit area disk. We study three cases of interest: 
an arbitrarily chosen narrow transmission band; the case of 
power limited nodes which transmit in disjoint narrow bands; 
and the case of nodes with high power capabilities that use 
of a wide transmission band. The choice of transmission band 
in the last two cases depends on the transmission distance 
and the physical characteristics of the channel, and is made in 
accordance with the waterfilling principle. 

We have shown that the amount of information that can be 
exchanged by each source-destination pair in an underwater 
acoustic network goes to zero as the number of nodes n goes 
to infinity, at least at rate n ~ 1 / a e~ w o ( -° ( ~ n ~ 1/a ^ . This rule 
is valid for the different scenarios in general, requiring only 



changes in the scaling constants. The throughput per source- 
destination pair has two different regions. For small n, the 
throughput decreases very slowly as n increases. For large 
n, it decreases as n~ l / a . Considering that 1 < a < 2 in 
an underwater channel, the available throughput for large n 
decays more rapidly than in typical radio wireless networks. 
However, typical node densities in underwater correspond to 
the small n regime. In a narrow band example with values of 
a(f) characteristic of an underwater channel, we showed that 
the upper bound on the throughput remains almost constant 
for densities of less than 100 nodes per km 2 . Most under- 
water networks have node densities in this range due to the 
applications for which they are deployed. 

Finally, we have pointed out some important characteristics 
of the underwater acoustic channel useful in future studies. 
For example, we could allow cooperation between nodes a la 
Ozgur et al [3] taking advantage of the distance-band separa- 
tion. That is, instead of performing time division between long 
and short transmissions, we could simply transmit in different 
bands that do not interfere with one another. This is important 
because acoustic transmissions have long propagation delays 
due to the speed of sound underwater (^1500 m/s), which 
reduces the usefulness of a time-division scheme. 
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